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ABSTRACT 

We derive the Hamiltonian describing Pauli-Fierz massive gravitons on a flat Friedmann- 
Robertson- Walker (FRW) cosmology in a particular, non-generic effective field theory. The 
cosmological evolution is driven by a scalar field $ with an arbitrary potential V ($) . The 
model contains two coupled scalar modes, corresponding to the fluctuations of $ and to 
the propagating scalar component of the Pauli-Fierz graviton. In order to preserve the full 
gauge invariance of the massless version of the theory, both modes have to be taken into 
account. We canonically normalize the Hamiltonian and generalize the Higuchi bound to 
FRW backgrounds. We discuss how this bound can set limits on the value of the Pauli- 
Fierz mass parameter. We also observe that on a generic FRW background the speed of 
propagation of the scalar mode of the graviton is always smaller than the speed of light. 
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1 Introduction 



The discovery of the recent cosmic acceleration in the late '90s has revamped the interest on 
theories that attempt to modify gravity at large distances. An interesting area of studies is 
to allow for a small, but finite graviton mass. By assuming the mass to be sufficiently small, 
the gravitational interaction should depart from the predictions of the theory of General 
Relativity only at large distances, or equivalently late times, that is at scales comparable 
with the graviton Compton wave-length. 

It was shown by van Dam and Veltman [1] and by Zakharov [2j that this is in fact not 
the case on flat Minkowski background. For instance, for a mass of the order of the current 
value of the Hubble parameter, modiflcations are already manifest at Solar System size, 
where the contribution coming from the mass would have been expected to be irrelevant. 
The discontinuity, known as vDVZ discontinuity, is only present for a flat background metric, 
for it was shown in [3l H] that the zero-mass limit is smooth when a non-zero cosmological 
constant is present. The removal of the vDVZ discontinuity in (Anti)-de Sitter space is due 
to the existence of two mass scales - the cosmological constant and the graviton mass 
m?. However, for positive values of the cosmological constant, a new pathological regime 
appears: if the mass of the graviton does not obey the inequality > 2H^, known as the 
Higuchi bound O |6] , the theory develops a ghost-like instability. 

In the light of the recent interest in potentially phenomenological implications of theories 
in which gravity is modifled in the infrared regime, it is interesting to discuss how such a 
bound can be extended to more realistic backgrounds than pure de Sitter space. It is therefore 
the goal of the present work to investigate the effects of a Pauli-Fierz mass term for the 
graviton when the background space-time is described by a flat Friedmann- Robertson- Walker 
metric (i.e., when the mass scale H is time dependent) and, in particular, to generalize the 
Higuchi bound to a FRW UniverseE 

In order to consistently account for the dynamics of the energy density that drives the 
expansion of the Universe, we will consider a Pauli-Fierz graviton coupled to a scalar fleld $ 
with a generic potential V ($). The background metric is characterized by a time-dependent 
Hubble scale H{t). In order to maintain the full coordinate reparametrization invariance of 
the massless theory, we have to consistently take into account the fluctuations of $ along 
with those of the graviton. This leads, in the scalar sector of the theory, to a complicated 
coupled dynamics. By following closely the analysis of [B], we flnd the generalized Higuchi 
bound to now depend - through the behavior of $ (t) - not only on the Hubble scale H{t) 
but also on its flrst time derivative H{t). This bound has to be satisfled at any time during 
the cosmological evolution for the theory not to develop any instabilities. If we assume the 
bound to be universal, that is to be independent on the species that drives the cosmological 
expansion, our result strongly constrains the value of the Pauli-Fierz mass. 

Furthermore we will show that a physical degree of freedom of this model displays a 
Lorentz-violating propagation, that however never develops into superluminality. 



*It is worth to mention that in a recent paper [7] a mechanism was described to modify the Higuchi bound 
in a de Sitter Universe by adding a suitable coupled ghost; in our work we will consider only canonically 
normalized scalar fields. See also |8] for the study of the behavior of a different model of infrared-modified 
gravity on FRW background. 
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The paper is organized as follows: in section 2 we will illustrate the Hamiltonian of the 
system describing gravity with a Pauli-Fierz mass term and a scalar field. In section 3 we 
will diagonalize the part of the Hamiltonian that describes the propagating tensor and vec- 
tor modes. In section 4 we will discuss the conditions under which the scalar sector can 
be brought to canonical form (details can be found in the appendices) and, as usual, the 
most interesting dynamics occurs there. The implications of this non-trivial dynamics will 
be discussed in section 5. 



2 Quadratic Lagrangian of Massive Gravity on FRW 

Our system consists of a scalar field $ canonically coupled to gravity, which we allow to 
slightly depart from standard GR. 

S = Sq + Sm ■ (1) 

5*0 is given by 

So^-Jd'x^/^ , (2) 

and Sm is the most general quadratic self-interaction of the metric that does not exhibit 
(ghost-like) instabilities and preserves the Lorentz invariance of the Einstein-Hilbert action. 
The term is known in literature as the Pauli-Fierz mass term: 

^ J d'xLm = y d^x^/^h^, [r^r" - r^gn v ^ (3) 

where cjf^u is the background metric that solves the equations of motion at zeroth order in 
the fluctuations h^u = g^y — g^u over said background. By construction the Pauli-Fierz term 
does not contribute to the equations of motion for the background since it is quadratic in 
the fluctuations, therefore the metric g^u is found by solving the standard Einstein equations 
of motion in presence of a scalar field 

Before proceeding any further it is worth to make few comments. Firstly, in the present 
paper we will study the properties of the model ([1]) in the spirit of effective field theory: we 
will be assuming that a UV completion does exist and reduces at low energies to a massive 
theory of gravity described by ([T]) with mass given by 

One might wonder why we constrain ourselves to a Lorentz-invariant mass term, while 
the background we are interested in is ultimately Lorentz-violating. As we will see in sec- 
tion [Sj the Lorentz-violating background will give rise to Lorentz-violating phenomena: in 
the high frequency regime, the speed of propagation of (the scalar component of) the grav- 
itational modes is smaller than the speed of light. Hence it might seem natural to allow 
for a more general mass term, one that will only preserve the symmetries of a FRW back- 
ground. The reason for not doing so is twofold: on one side, we would like to make contact 
with the literature on PF massive gravity - the first and foremost motivation of this work 
was to extend the Higuchi bound to cosmological backgrounds. On the other end, we want 
to analyze what kind of Lorentz-violating effects arise because of a background that only 
preserves a subgroup of the symmetries of the theory; considering a more general quadratic 
self-interaction would only entangle those effects coming from a non-maximally symmetric 
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background to those generated by a Lorentz-violating modification of tlie Einstein-Hilbert 
action. 

Tlie action can be rewritten in tlie ADM formalism as: 



So = J d^x {^n'^gij + H $ + [iV £° + A^. r ] } 



(4) 



where £^ = 11^- r° and 



TZ' - V and 
V9R + 



y/9 



1 



m 



T = Ild'^. (5) 
The following quantities, known as lapse and shift function, are usually not dynamical 



N = ll^-g^\ N, = go,. (6) 
The canonical momenta, and 11, are related to the canonical coordinates, gij and $, via 



2N 



9ij 



^/9 

NU 
\f9 



(7) 



The PF mass term can be rewritten in this formalism as 

,2 



J d'^x a [hij hij — ha hjj — 2 a? NiNi — Aa? n hii\ 



The background solution for $ is denoted as <l>(t) and depends only on time since we 
are interested in a cosmological background, which by construction does not depend on 
(jSj). This has the form of a flat Friedmann- Robert son- Walker, that is dsg = (^^j^dx'^dx^ = 
-dt^ + it) dx^. 

We then consider small perturbations over this background. We will study the dynamics 
of perturbations over this background. It is known that in theories of massive gravity there 
is an unexpected very low strong coupling scale, usually of order A5 = (Mpm^)^/^ and 
above such scale perturbative theory can not be trusted; the following analysis is therefore 
considered only within the range of validity of the effective theory, that is for wavelengths 
smaller than the Hubble radius and longer than the strong coupling scale. The study of the 
nonlinear effects on top of this theory is part of a future follow-up of this work. 



I^We will conservatively estimate in the present work A5 to be the strong coupling scale: the precise 
evaluation of this scale will be the subject of a future work. 
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Following the notation of [6] the fields can be written as 

Tiv = T^ij + pij = -2 Ha 5'^ + p'^ 

N = 1+n, $ = l> + (^, n = n + 7r. (9) 

Ni is already a first order quantity, since its background value vanishes. 
At quadratic level, the Lagrangian has the following form 

LQ = p'^ hij + 7T^ + nSl + S'Q + Ni 4 + Lm (10) 

where the subscripts L and Q denote linear and quadratic order of the quantities under 
consideration. 

In the absence of the graviton mass term, both n and Ni appear linearly in the Lagrangian. 
Not being propagating degrees of freedom, their equations of motion enforce four constraints 
on the physical propagating modes. Hence the total number of physical degrees of freedom 
is two for the graviton, as expected from the Lorentz representation of a massless tensor 
field, and one for the scalar field. 

When the mass term is introduced, the dynamics changes. The absence of time derivatives 
for Ni still suggests that this mode is not propagating, but the mass term introduces a 
quadratic term for A^j: its equation of motion is therefore an algebraic constraint for Ni 
itself. Two additional modes begin to propagate in this model, a transverse vector and 
a scalar field, accounting to five degrees of freedom for the graviton as demanded by the 
Lorentz representation of a massive tensor field. 

The Lagrangian reduces to the following, once Ni is integrated out 

2 / ipi \ 2 2 

L = p'^ hij + 7i^ + n {£l + — h,) +£l- - ^ {Kj hij - ha hj^) . (11) 

Following [6], the expression can be furthermore simplified through the redefinitions: 
hij — )■ a^/'^hij, p^^ — )■ a~^/^p*-', n — a'^/'^n, — )■ a"^/^ vr — )■ a^/^vr. We also define 
A = didi/a^ (note that A depends on time) and rescale 11 — )■ a^fl (z.e., fl = $). 

The background equations are 

a2 12 ^ 6 ' 
4 ' 

11 = -?,HI{-V' (12) 

where V = V{^) etc. 

Finally, we exploit the 5*0(3) invariance of the background to decompose hij in tensor, 
vector and scalar modes: 

5. ^a„,5,^i ,.,||,. (13) 

and analogously p'^K Through this definition we can now study the scalar, vector and tensor 
components of the action separately. 
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3 Tensor and Vector Modes 

The Lagrangian for the tensor modes hj^ takes the form 



^tensor 



+ - h^' 
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4 / 



(14) 



that is easily brought to canonical form by defining p^* — )■ pij/\/2 + 5H qij/{2\^) and 
hJJ^ — )■ \/2qij. The final form of the Lagrangian is thus 



-^tensor Pij Qij 



9 3 • 



(15) 



and we see that tensor modes have always well-behaved kinetic terms and, as usual, obey 
the same equation of motion of minimally coupled scalars. 
The Lagrangian of the vector modes, hj, can be written as 



^vector 



2p, ( 1- Pi + lli (m^ + 4H'-16H'^ + n^ ] q, + 





A \ 1 


Hpi (^-5 + 


8 — Ui 







(16) 



where we have defined pi = \/—didi p\ and q, = yj—didi h\. As in [B], we perform the following 
canonical transformation 



Pi 



AH Pi - (m^ -6H^) qi 
2m 



3Hqi + 2pi 
2 m 



(17) 



that brings the Hamiltonian in canonical form 

_ . ri 1 

-^vector — PiQi ~ '^PiPi ~^ 



-A + m'-'^-H^ + ^-H] q. 



Like for the tensor modes, also the vectors always present a well-behaved Hamiltonian (as 
long as m 7^ - of course, for m = the vector modes turn into purely gauge modes). 



4 Scalar Modes 

Our starting Lagrangian has the form 

Lscalar = P^ + + TT - {p\ p\ TT, h\ h\ (p, Tl) . (19) 

Since n appears in the Lagrangian as a Lagrange multiplier, its contribution disappears 
once the corresponding constraint equation + m? a~^/^ = is integrated. For scalar 
modes, such an equation reduces to 

i^^ {h^ + h^) - A - 2 H {p^ + p^) - V' if - Un = , (20) 
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where we have defined the quantity 



z/2 = - 2 + n72 = - 2 - 2 H = - 2- , (21) 

a 

that reduces to the Deser and Waldron parameter z/^ in the hmit if — )■ [6]. 

We use eq. ( |2T1) to ehminate p* from our Lagrangian. Through several integrations by 
parts the Lagrangian can be written as 

Lscalar = Po % + Pi Ql - ^-1 {Po, Pi, ?0, ^1, (22) 

where we have defined the new canonical variables 

qo^h'-j, q,^v-^hr (23) 

The variable go is analogous to the one defined in [6]. The variable gi corresponds (modulo 
an overall factor a (t)) to the Mukhanov variable v [10], that is the canonically normalized 
scalar degree of freedom in a Universe filled by a scalar field. These are the two fundamental 
scalar degrees of freedom of our system. 

The variable is not a dynamical degree of freedom, as it lacks its own canonical 
momentum p*. Variation of the action with respect to /i* gives an algebraic equation for 
whose solution reads 



3m^ jy2 



-mV^ qQ + [Upi + V' qi) - AH A{po - H go - ^ qi] 



(24) 



Note that the above solution is singular for z/^ = 0. As discussed in [HI El [12], the critical 
line z/^ = corresponds to a partially massless theory, where /i*, appearing only linearly in 
the Lagrangian, is a Lagrange multiplier. As a consequence, the corresponding equation for 
z/^ = provides one additional constraint that removes the scalar mode of the graviton from 
the propagating degrees of freedom. The result is that - at quadratic level - the graviton 
is left with only four helicities (ib2, ±1). As we will discuss in section 5, the possibility of 
living on the critical line z/^ = is much less interesting in the generic FRW background 
considered here, since the line z/^ = can be crossed, during the evolution of the Universe, 
only when system is already in a unstable region of the parameter space. 

Let us now go ahead and consider a generic z/^ 7^ 0. By plugging the result back 
into (!22|) we obtain the final Hamiltonian. Since its form is rather complicated, here we will 
express it in the following compact notation 

Lscalar = P^ ■ q - (^-p^-K-p+]^q^-M-q+p^-V-(^ (25) 

where we have defined the vectors = {po, pi) and q^ = (go, gi) and where the explicit 
expression of the matrices K, M and V are given in Appendix A. 
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-2E n— go A + ' ' (26) 



By performing a series of canonical transformations (details in Appendix A) we can bring 
the Hamiltonian of our system to the form 

V - ^^P2 i/M^^ + 2A^) Pi QiAQ, 

- 24W^°"2 — QoAQo + ^ ^ + 

where we have defined a new function of the background quantities 

X'^ = m^-2H^, (27) 

and where the elements of Mjj can be obtained by setting A = in eqs. ( 1461H8|) in Appendix 
A and are complicated functions of the background quantities. 

Even if eq. (|26ll does not yet correspond to the Hamiltonian in its canonical form P^/2 — 
QAQ/2 + V{Q), it already shows an important peculiarity of our system. We see that the 
sign of the coefficient of the kinetic term of the gravitational mode Pq that of z/^. This 
observation allows us to obtain a first extension of the Higuchi condition > 2 to a 
generic cosmological background. In a cosmological background, we notice that a necessary 
condition for the positiveness of the Hamiltonian reads z/^ > 0, i.e., 

m^>2H^ + 2H. (28) 

If the above inequality were to be violated, then the kinetic term of the degree of freedom 
Qq would have the wrong sign, signaling an instability of the theory. 

From now on we will assume that the inequality ( 128|) is satisfied. Then, we can canonically 
normalize the variable Qo through the redefinition Pq — ^ -Po7 — iQo, Qo Qo/l where 
7 = 2 Hu^/S/ (mA^). We thus get the final result 

-n^goPi + ^^^^^ (29) 

where Mqo has a new form, given in eq. (H9|) in Appendix A, Mqi is given by the same Mqi 
of eq. (H7|) divided by 7, and Mn is unchanged. 

Eq. (l29l) is our final expression for the Hamiltonian governing this system. In principle, 
one can perform a canonical transformation to eliminate the QqPi- This is however not very 
illuminating, so we will relegate to Appendix B the discussion of the strategy that allows to 
eliminate such a term. 

Before discussing our results, we check that in the limit of de Sitter background {i.e., H, 
V' 0) we recover the result of |6]. Indeed, in this limit the Hamiltonian reduces to 

nf = ^ + ^ + ^ (-A + m'-^nA go + % (-a + v"~^hA q, (so) 



2 2V A J ^ 2 ^ 

Of course, this canonical form can be obtained only if > 0, that is when the Higuchi 
bound > 2 is satisfied for a Pauli-Fierz graviton on a de Sitter background. 
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5 Discussion 



The Hamiltonian for a massive graviton on a cosmological background governed by a scalar 
field has two remarkable properties. First, in order for the kinetic term of the gravitational 
modes to have the right sign, the system must obey the inequality > 0. Second, since 
the background breaks the isometry group down to 5*0(3), the scalar mode of the graviton 
turns out to inherit a Lorentz non-invariant dispersion relation. 

We will assume the condition z/^ > to be met in the following discussion. 

For wavelengths shorter than the characteristic timescale of the background and shorter 
than the Compton wavelength of the graviton, the mode Qq follows the equation 

Qo-Q + ^^) AQo = 0. (31) 

The coefficient of the Laplacian in the previous equation has to be positive to prevent the 
generation of exponentially increasing solutions. The explicit expression of this coefficient is 

1 2A2 AH m'^-2H^-2H/3 , , 

o + ^ = 1 + ^ = 2 ■ (32) 

Before discussing the stability of the model, let us note that the speed of propagation 



of these modes is y 1 + 4i^/3z/^. Since < 0, it is guaranteed that these modes do not 
experience superluminal propagation. Interestingly, [13] has observed that the speed of 
propagation of gravitational modes cannot be smaller than (at least) (1 — 2 x 10~^^) c. Were 
this bound violated, gravitational Cherenkov radiation would have depleted the population 
of high energy cosmic rays that we currently observe. Unfortunately, this bound cannot 
be directly applied to our Pauli-Fierz gravitons, in fact it relies on the tensorial coupling 
of the graviton, while we were interested in the couplings of the scalar component of the 
Pauli-Fierz graviton. Moreover, the Pauli-Fierz theory is known to be strongly coupled at a 
very low scale of the order of m [Mp /m)^^^ [l4j, while the bound mentioned above relies on 
the emission of gravitons with energies of the order of 10^^ GeV. It would be interesting to 
study whether a bound, analogous to that of [I3], could be applied to the scalar component 
of the Pauli-Fierz graviton. 

The requirement of positivity of eq. (132!) . i.e., > 2 + 2H/3, provides a more 
restrictive condition than i/^ > 0. By considering the equation of state parameter of the 
background w with —1 < w < 1, H and are related as H = —3 (1 -|- w) H'^/2. The 
condition z/^ > can be written as m? / H'^ > — (1 + 3w), whereas the condition (1/3 + 
2A^/3z/^) > is equal to m'^/H'^ > 1 — w. It is easy to see that the second condition is 
always stronger than the first one and that they coincide only on a de Sitter background. 
As a consequence, the lower limit on the mass of the graviton on a cosmological background 
is given by 

> (1 - w) . (33) 

We note that, since the requirement of positivity of eq. (l32l) is more restrictive than the 
bound z/^ > 0, the "partially massless" regime discussed in [HI [12] can not be attained as a 
consequence of the cosmological evolution. While for a de Sitter background is a constant, 
in our case it is time-dependent, so z/^ might cross zero during the cosmological evolution. 
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However, when this happens, the quantity ( !32l) is aheady negative, so the theory is aheady 
in an unstable regime hence such a hmit can not be trusted. 

Assuming that the constraint f l33p is independent of the kind of matter that is driving the 
cosmological expansion, we can infer a phenomenological bound on the mass of the Pauh- 
Fierz graviton. Indeed, we know that the Universe had a standard cosmological evolution 
at least from the time of nucleosynthesis, when it was in a thermal bath at a temperature 
of ~ 10 MeV. If we assume the mass of the graviton to be constant during the cosmological 
history of the Universe, then eq. ( l33l) constrains the Pauli-Fierz mass to be larger than the 
value of the Hubble parameter at nucleosynthesis, m > 3x 10~^^ eV. Such a mass corresponds 
to a Compton wavelength of the order of tens of thousand of kilometers, which is ruled out 
since the solar system dynamics is well described by newtonian gravity. Of course, the limit 
m > 3 X 10~^^ eV does not imply that a massless graviton is ruled out. Indeed, as we have 
discussed in section 2, we expect our linearized theory to break down for distances smaller 
than the scale at which the theory becomes strongly coupled. For values of m so small that 
the strong coupling scale is comparable to the Hubble radius (and a fortiori for a massless 
graviton), we will not be able to trust our linearized analysis. 

To conclude, we have derived the canonical Hamiltonian describing Pauli-Fierz massive 
gravitons on a FRW background driven by a scalar field. The scalar sector contains a mode 
with Lorentz non-invariant dispersion relation. The requirement that the momentum part of 
the Hamiltonian is positive definite induces the bound ( 128|) . whereas the (more restrictive) 
bound (133|) prevents the development of instabilities through the generation of large gradients 
of the fields. Our formulation provides a setting for the study of the presence or the absence 
of the vDVZ discontinuity, as well as the existence of strongly coupled regimes, on a FRW 
background. We hope to go back to these points in a forthcoming publication. 

Note added 

While we were at the final stages of the writing of this paper, the work [TS], whose 
subject partially overlaps with ours, was posted on the archive. The generalized Higuchi 
bound obtained with the Pauli-Fierz choice of parameters in [15] does not agree with our 
eqs. (128|1 or ( |33|) . due to different nature of the starting Lagrangians. 
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A Canonical Transformations 

In this appendix, we show explicitly the canonical transformations that allow to bring the 
Lagrangian fl2Sl) to the Hamiltonian (^E^. Our procedure follows rather closely the one of [H]. 
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Our starting point is the Lagrangian (125 



where 



i^scaiar = ■ q - ( ^ ■ K ■ p + ^ ■ M ■ q + ■ V ■ q ] , (34) 



q _ 4A _|_ JL f'i — 

^ — \ n /q _ 2A\ 1 I J l-J^J 



n 



3 -H 

--2 ^ 
^ 3;/2 



12 /f2 



An 



+ 



v 

2H 



_ AV 



12 /f2 1,2 



_ 2 _ 

3 m 

n2 A 

6// 1/2 



A2n 



(36) 



M 



^2,tt2 4AH^ 4^2j^ 2 a , rn? v"^ ( 2A m'' \ i7/ i / 4 g A 4 g A \ rr 

T4-7rE^)v"+(4^-^-^)n -A + y"- (1-^) + ^:^ (£122^ + AH 



(37) 



First, we rewrite it in the form 

1 , . . . . 1 



i^scaiar = p' ■ q - [ - {p - h ■ q)' ■ K ■ {p - h ■ q) + - q' ■ M ■ q + {p - h ■ q)' ■ V ■ q ] (3^ 



where the matrix 



H n/2 

12 

2A^ 



n/2 -V'/U-3H + 



(39) 



is chosen in such a way that the matrices M = M + h'^ ■K-h+h'^ ■V + V'^ -h and V = V + K-h 
acquire a simpler form. Note that we have defined the new time-dependent function 

X'^ = m^-2H^. (40) 

The explicit expression of the matrices M and V is such that 





M + h 



121/2 



-2 3 rr I A A2 fr 4A-3n2 



■m _i_ 3 rr I A2 TT 

'2H 2 y 3 1^2 ^ 4^2 

m2A2n _ /J [^3 _|_ 2n2 



y ^ _ _ 

' " 12 H2 1/2 If \2"^^ 

We perform the following canonical transformation 

p = p + h- Q + h- a- P , q = Q + a-P , 

[4A4 (3m2-2 A-9H2)-(l2y'+54J? ft) HnA2-3 (m2+2_H"2^ n*] fv / ^2 
a= I Y _m2 A6 ^V-^ 

n/A2 
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followed by the simple canonical transformation Pq ^ Qo, Qo ^ —Po- The Lagrangian now 
has the same form as in eq. (1251) where now 

V 1 

V = — Tin? — _+ — — _ _ u _ m 
' -2HnuyX^ -r/n-/7 (3/2 + 2 nVA2) ' ' 

Moo Moi 
Moi Mu 



M 



with 

Moo = _ 4gM^y2A^) ^^ 3H' ^^q x^u' V" - 32 V"+ (42) 

- lei/A^n (i2A2z/2 + n^) v' + s (4m2 -gi^^^ a^ + 12 (m^ - 40 i/^) A^ n^] , 
3 n 

Mil = -A + ^ [Si/A^ + 2i7n) +m^n + 6i/2n^] , (43) 

Moi = [2 {V' + HU) X' + SHUu'] , (44) 

We can eliminate the diagonal terms of the matrix VhjP^P + A- Q, Q ^ Q where 
the matrix A is defined as diag(— Voo/-f^oo; —Vu/Kn). This leaves the matrix K unchanged, 
while V simplifies to 

^=(-2HLyx^ o) ' ^4^) 



and the matrix M has components 

{jy^ + 2X^^ 
m? A^ ' IV? X^ 



(z/2 + 2A2) HU r 9 9 9 / 9 9M 

Moo = -4 ^-A + Q^-[2m'u'-X' {m' + 3H^)]V' + 



U (4m'-9H') A^ + 2 (llm^-39H') A^ff+ 
+ (19 - 44 H'^ + 12 H^) + 6 II^] , (46) 

Moi = Mio = -2^ (5Z/2-4A2) (6A^ + 35A2ff + 5n^) , (47) 

^„ 9 , 2/71^-n 3 (m^-10ff^)ff 1 (m^ + /7^) ^ 
Mn = -A + y --i/ -^^^+8 2 A^ " ^^^^ 

As discussed in section 4, if z/^ > then we can bring the P^ terms to canonical form 
by redefining Po — ^ -Po7 — iQo, Qo Qo/l with 7 = 2 Hu\/3/X^. The Hamiltonian then 
keeps the same form as above, however the coefficient of the Pq term is now just 1/2 {i.e., 
Koo = 1), Moo 

[8 (4m2 - 9H^) X^ + 12 (m^ - AOH^) X^ 



32 X^ 



- 6 (m^ + 47/72) A^n^ - (I7m^ + 20m^ - 108H^) ^-AH^ W] 
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whereas Mqi goes to M01/7 and Vio goes to Vio/7- 



B Removing the Qo Pi term 

In this appendix we show how to formally bring the Hamiltonian derived in appendix A to 
the canonical form ^p^ ■ p + ■ Mgn ■ <?• We start from the expression 

'H = ^pT .p^^qT . M .q + p^ -V -q , V = 

where v{t) is an arbitrary function of time and M is an arbitrary time-dependent 2x2 
"mass" matrix. We look for a transformation of the form 

p^A-p + B-q, q^A-q. (51) 

that can be shown to be canonical if the matrices A and B are such that ■ B is symmetric. 
In order to leave the kinetic term untouched, the matrix A has to be orthogonal. Let us 
write it as 

^_ ( cosa(t) sina(t) 
\^ — sina(t) cosa(t) 

The matrix B has to be equal to A — V ■ A to eliminate the cross term p^ -V ■ q from the 
Hamiltonian. The function a(t), finally, is determined by the requirement that A^ ■ B is 
symmetric. It is easy to prove that such condition is verified if the function a{t) satisfies the 
differential equation a{t) = —v{t)/2. Both A and B are then determined up to an integration 
constant. By using the equation that relates A to 5 we find that the Hamiltonian can be 
brought to the final expression: ^p^ ■ p + ■ Mgn ■ q with Mgn = MA — A'^ V'^ V A — 
A^A-A^VA. 




References 

[1] H. van Dam and M. J. G. Veltman, "Massive and massless Yang-Mills and 
gravitational fields," Nud. Phys. B22 (1970) 397-411. 

[2] V. I. Zakharov, "Linearized gravitation theory and the graviton mass," JETP Lett. 12 
(1970) 312. 

[3] M. Porrati, "No van Dam-Velt man- Zakharov discontinuity in AdS space," Phys. Lett. 
B498 (2001) 92-96, h ep-th/0011152, 

[4] I. I. Kogan, S. Mouslopoulos, and A. Papazoglou, "The m limit for massive 
graviton in dS(4) and AdS (4): How to circumvent the van Dam-Veltman-Zakharov 
discontinuity," Phys. Lett. B503 (2001) 173-180, h ep-th/0011138, 

[5] A. Higuchi, "Forbidden mass range for spin-2 field theory in de Sitter space-time," 
Nucl. Phys. B282 (1987) 397. 



12 



S. Deser and A. Waldron, "Stability of massive cosmological gravitons," Phys. Lett. 
B508 (2001) 347-353, |hep-th/0103255| 

G. Gabadadze, A. Iglesias, and Y. Shang, "General Massive Spin-2 on de Sitter 
Background," 10809 . 29961 

C. Deffayet, "Cosmological perturbations of brane-induced gravity and the vDVZ 
discontinuity on FLRW space-times," Phys. Rev. D71 (2005) 103501, g r-qc/0412114] 

R. L. Arnowitt, S. Deser, and C. W. Misner, "The dynamics of general relativity," 
gr-qc/0405109, 



V. F. Mukhanov, H. A. Feldman, and R. H. Brandenberger, "Theory of cosmological 
perturbations. Part 1. Classical perturbations. Part 2. Quantum theory of 
perturbations. Part 3. Extensions," Phys. Kept. 215 (1992) 203-333. 

S. Deser and R. I. Nepomechie, "Gauge invariance versus masslessness in de Sitter 
space," Ann. Phys. 154 (1984) 396. 

S. Deser and A. Waldron, "Partial masslessness of higher spins in (A)dS," Nucl. Phys. 



B607 (2001) 577-604, [h^^h/0103198[ 



G. D. Moore and A. E. Nelson, "Lower bound on the propagation speed of gravity 
from gravitational Cherenkov radiation," JHEP 09 (2001) 023, |hep-ph/0106220: 



N. Arkani-Hamed, H. Georgi, and M. D. Schwartz, "Effective field theory for massive 
gravitons and gravity in theory space," Ann. Phys. 305 (2003) 96-118, 



hep-th/0210184| 



D. Bias, D. Comelli, F. Nesti, and L. Pilo, "Lorentz Breaking Massive Gravity in 
Curved Space," 10905.16991 



13 



